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2 I Abstract 

in ■ 

O I The charged particle multiplicity in hadronic three-jet events from Z decays is 

><; ■ investigated. The topology dependence of the event multiplicity is found to be 

^ . well described by a modified leading logarithmic prediction. A parameter fit 

Q-i. of the prediction to the data yields a measurement of the colour factor ratio 

^ I Ca/Cf with the result 

.> ; Ca/Cf = 2.261 ± 0.014stat. ± 0.036exp. ± 0.066theo. 

\-{ I in agreement with the SU(3) expectation of QCD. The quark- related contribu- 

tion to the event multiplicity is subtracted from the three-jet event multiplicity 
resulting in a measurement of the multiplicity of two-gluon colour-singlet states 
over a wide energy range. The ratios r = Ngg{s)/Nqq{s) of the gluon and quark 
multiplicities and r^^^ = Ngg{s) / Ng^{s) of their derivatives are compared with 
perturbative calculations. While a good agreement between calculations and 
data is observed for r^^\ larger deviations are found for r indicating that non- 
perturbative effects are more important for r than for r^^\ 
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1 Introduction 

Abundant particle multiplicity is one of the most obvious properties of hadronic final 
states in e'^e~ annihilation. The reason for the large multiplicity of particles produced 
in hadronic events is directly rooted to one of the fundamental properties of the strong 
interaction: the confinement of quarks and gluons into hadrons. Unlike leptons, quarks 
produced in a high energy reaction are not able just to separate without further inter- 
actions in which the colour charges are balanced and free colour singlets are formed. 
Ironically, it is the very same property of QCD which makes it impossible to predict the 
average number of hadrons produced in such an event. However, the hypothesis of Local 
Parton-Hadron Duahty (LPHD) allows the assumption to be made that the hadronic mul- 
tiplicity of an event is related only via a normalisation factor to the partonic multiplicity 
at a given virtuality cut-off Qq, which then is a perturbatively calculable quantity. 

The multiplicity of qq colour-singlet systems is well understood in this way. Theoretical 
predictions of the energy dependence of the multiplicity produced in such systems [1- 
3] have been confirmed by experimental data to high accuracy [4-6]. It is possible to 
predict the multiplicity of two-gluon colour-singlet systems perturbatively in a similar 
fashion, but experimental verification of this quantity is scarce. Two-gluon systems are 
experimentally difficult to produce and are observed so far only in decays of resonances 
with relatively low mass. The ratio of the multiplicity of quark and gluon systems, r, 
has been subject to theoretical studies presented in Sect. 3. It is expected that r at high 
energies resembles asymptotically the colour factor ratio Ca/Cf- Cp and Ca, which are 
the Casimir eigenvalues of the triplet and octet representation of SU(3), can be interpreted 
as effective colour charges of the triplet quarks and the octet gluons. However, r is known 
to have large corrections in higher orders of the perturbative expansion and moreover, is 
supposed to be affected by non-perturbative effects. Therefore, the ratio of the derivatives 
of the quark and gluon multiplicities with respect to energy, r^^\ has been suggested in [7] 
as a better suited observable. 

The next slightly more complex class of events which can be studied in e'^e~- 
annihilation are events with three jets, which occur if in a qq-event a gluon is radiated 
with a sufficiently large transverse momentum. Three-jet events are unique in the sense 
that they provide quark jets as well as gluon jets and therefore provide the opportunity 
to study the properties of gluon fragmentation once the gg-contributions are understood. 
There have been numerous experimental studies where the properties of identified quark 
and gluon jets are compared (e.g. [8-12]). A crucial question arising in this context is 
how the environment of a three-jet event affects the properties of a quark or gluon jet 
compared to the unrestricted jets of a qq- or two-gluon event in terms of phase space 
restriction or coherent radiation. Studies have been made to verify experimentally sev- 
eral proper energy scale like variables which account for these differences by applying 
the appropriate evaluation scales for the jets in three-jet events (e.g. in [9,11]). Coher- 
ence effects in the particle production in three-jet events have been subject to a recent 
study [13]. In parallel with the publication of an analysis of the charged multiphcity of 
symmetric three-jet events applying a rather phenomenological approach [14] , a theoret- 
ical prediction of the multiplicity of three-jet events has been published [15,16], which is 
derived in the colour-dipole picture of the Modified Leading Logarithmic Approximation 
(MLLA) and takes into account the effects of colour coherence and phase-space restriction 
in a stringent way. This prediction has been applied to study the topology-dependence 
of the event multiplicity of symmetric three-jet events in [17] and later in [18]. 
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In this paper this prediction is not only apphed to the charged multiphcity of three-jet 
events of symmetric topologies, but events with more general topologies are also con- 
sidered. In the following section the experiment, data selection and the classification of 
three-jet event topologies is laid out. Sect. 3 deals with the relevant theoretical predic- 
tions, before in Sect. 4 the measured multiplicities are discussed. The following section 
deals with the preparation of the predictions which are then compared to the data in 
Sect. 6. A parameter fit to determine the colour-factor ratio Ca/Cf is performed and 
the systematic uncertainties of this measurement are discussed. In Sect. 7 the prediction 
is used to extract the multiplicity of two-gluon colour singlet systems from the three-jet 
event multiplicity by subtracting the quark contribution. With the measurements of the 
charged two-gluon multiplicity Ngg the ratios r and r^^^ are evaluated and the possibil- 
ity of measuring the ratio of the second derivatives, r^'^\ is studied. The summary and 
conclusions are presented in Sect. 8. 



2 Data and data analysis 



variable 


cut 


P 


> 0.4 GeV 


^polar 


20° - 160° 




< 5 cm 




< 10 cm 


-^track 


> 30 cm 




< 100% 


-Ehpc 


0.5 GeV - 50 GeV 


-E'emf 


0.5 GeV - 50 GeV 


-E'hac 


1 GeV - 50 GeV 



Table 1: Selection cuts applied to 
charged-particle tracks and to calorimeter 
clusters. 



variable 


cut 


general events 


phemisph. 
charged 


> 0.03 ■ 


ptotal 
charged 


> 0.12- 


-^charged 


> 5 


^sphericity 


30° - 150° 




45 GeV 


three-jet events 




> 355° 


-E'visible/jet 


> 5 GeV 


^charged /jet 


> 2 




30° - 150° 



Table 2: Selection cuts applied to gen- 
eral events and to three-jet events. 



In this paper the hadronic events from Z decays recorded by the DELPHI experi- 
ment in the years 1992-1995 are analysed. The Delphi detector was a hermetic collider 
detector with a solenoidal magnetic field, extensive tracking facilities including a micro- 
vertex detector, electromagnetic and hadronic calorimetry as well as extended particle 
identification capabilities. The detector and its performance are described in detail else- 
where [19,20]. 

In order to select well-measured particles originating from the interaction point, the 
cuts shown in Tab. 1 were applied to the measured tracks and electromagnetic or hadronic 
calorimeter clusters. Here p and E denote the particle's momentum and energy, i?poiar 
denotes the polar angle with respect to the beam, is the distance of closest approach 
to the interaction point in the plane perpendicular to {xy) or along [z) the beam, respec- 
tively, I/track is the measured track length. -Ehpc (-Eemf) denotes the energy of a cluster 
as measured with the barrel (forward) electromagnetic calorimeter, and E'hac the cluster 
energy measured by the hadronic calorimeter. 

The general event cuts shown in Tab. 2 select hadronic decays of the Z and suppress 
background from leptonic Z decays, 77 interactions or beam-gas interactions. Further 
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reduction of background to a negligible level is achieved by the jct-sclection cuts given 
also in Tab. 2 . The cut variables are the visible charged energy, -^charged -^charge'? ' 
observed in the event or in each event hemisphere, respectively. Event hemispheres are 
defined by the plane perpendicular to the sphericity axis. The polar angle of this axis 
with respect to the beam is "(^sphericity and Ai'chargcd is the observed charged multiplicity. 
Events are discarded, if they contain charged particles with momenta apparently above 
the kinematic limit. 

In addition to the event selection above, a procedure to tag events with initial 6-quarks 
is applied. Details of the tagging procedure, which is based on the combined information 
of several observables, can be found in [21]. A cut on the tagging variable A by demanding 
A < 1.5 is used to reject 6-events leading to a hght quark event purity of ~ 92%. Both, 
the complete set of accepted hadronic (udscb-) events as well as the set of anti-tagged 
light (udsc-) quark events are analysed for cross-check reasons. 

In the accepted events three jets arc then reconstructed using the angular ordered 
Durham algorithm [22] taking into acount reconstructed momenta of charged particles 
and of neutral particles seen in the electromagnetic calorimeter. No cut-off variable for the 
cluster algorithm is used, every event is clustered into three jets. For cross-check reasons, 
the Durham [23], Cambridge [22] and Lucius [24] jet clustering algorithms have been used 
alternatively^. As three-jet events have to be planar due to momentum conservation, the 
three reconstructed jets are projected into the event-plane, which is defined by the first 
two eigenvectors of the sphericity tensor [25] . 

The three-jet event quality requirements, also shown in Tab. 2, assure well-measured 
jets. Here 9i denotes the angle between the two jets opposite to jet i, 'djet the polar angle 
of a jet and -Evisibie/jet the total visible energy per jet, taking into account all tracks 
and information from the electromagnetic calorimeter assigned to this jet. Assuming 
massless jet kinematics the topology of an event can be characterised by the three angles 
^1,2,3 between the jets. The inter-jet angles are ordered according to their size with 6i 
being the smallest and 63 the largest angle. This reflects the wide-spread convention due 
to which jets are numbered according to their energy with jet 1 being the most energetic, 
while the inter-jet angles are numbered according to their opposing jet. With massless 
jets this implies the numbering of the angles given above. 

The intcr-jct angles of the projected jets add up to 360°. Therefore giving two of the 
three angles is sufficient to fully determine the topology of an event. Here, 62 and 6*3 are 
used to define the binning. The angular intervals used to define bins of different topology 
classes are listed in Tab. 3. The allowed values for 9^ depend on the value of 63, as ^2 < ^3 
and 62 > 180° — ^^3/2, the latter being a consequence of 62 > Oi. 

Besides these general classes of topologies, classes of symmetric topologies are selected, 
where symmetric events are defined by the requirement of the two larger inter-jet angles 
of an event, 62 and 63 to be equal within small tolerances, i.e. 9^ < 62 + A9. An 
angular tolerance less than 5° is not meaningful, as Monte Carlo studies showed that the 
differences between the inter-jet angles in the partonic and in the hadronic state scatter 
with a standard deviation of ~ 5°. Therefore = 5° is chosen here. Since 63 has to be 
larger than 120°, the binning is completely determined by the choice of A^. Symmetric 
events have been used in several previous studies (e.g. [8-11,14,17,18]), exploiting the 
special property of these events. With jet 2 and jet 3 being a quark and a gluon jet 
with exactly the same topological environment, a similar energy scale is present. Here 
symmetric topologies, which are also contained as a subset in the general event topologies, 

^It is in the nature of the Cambridge algorithm, that depending on the structure of an event not every number of jets can 
be resolved. The small number of events where three jets could not be resolved are discarded when using the Cambridge 
algorithm. 
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9 


2 




ds 


98° 


- 100° 


119° - 


123° 


120° — 


130° 


100 


— 102 


126 — 


"1 OOO 

128 


130° - 


140° 


102° 


- 104° 


128° - 


133° 


140° - 


145° 


104° 


- 106° 


133° - 


138° 


145° - 


150° 


106° 


- 108° 


138° - 


143° 


150° - 


155° 


108° 


- 111° 


143° - 


148° 


155° - 


160° 


111° 


- 115° 


148° - 


155° 


160° - 


165° 


115° 


- 119° 











Table 3: The bins in 62 and 63 

are studied also separately. Due to the additional constraint on symmetric topologies, 
giving only one inter-jet angle is sufficient to fully describe the topology of a symmetric 
event. For consistency with previous studies 9i is used for this purpose. The reduction to 
only one independent variable considerably simplifies most correction procedures. Note 
that the symmetric events are not fully contained in the general event topologies listed in 
Tab. 3, as the minimum value for 61 allowed by the binning in Tab. 3 is 40°. In symmetric 
events jet 3 is produced with maximum energy, therefore smaller opening angles 61 can be 
accepted. However, Monte Carlo studies showed that the resolved hadronic jet structure 
does not refiect the partonic structure of an event, if 9i is smaller than ~ 20°, which 
therefore gives a lower boundary for 61. 

For each {62, ^3)-bin of the general topologies, and each ^^i-bin of the symmetric event 
sample, the multiplicity distribution is measured. To correct these distributions for de- 
tector acceptance, a matrix correction is used. For each angular bin a matrix Mmn is 
calculated by tracing in a Monte Carlo simulation how many particles m in an accepted 
event with n detected particles have been generated. In order to obtain the generated 
number of charged particles, particles with lifetimes shorter than 10~^s, like and A, 
have been forced to decay. is then applied to the measured distributions. Since gen- 
erated events which fail to fulfil the selection criteria are not considered when calculating 
Mmn and these events are biased towards low multiplicities, the correction is slightly 
overestimated. In order to correct for this overestimation, the multiplicity distributions 
are multiplied by the ratio of the multiplicity distributions of all generated events with 
the multiplicity distributions of all generated and not rejected events for the respective 
angular bin. While the application of increases the mean of the distributions by 
~ 30%, the multiplicative second correction results in a reduction by only ~ 4%. 

When correcting the multiplicity distributions of udsc-eyents, the matrix corrected 
multiplicity distributions still contain some contributions due to mistagged 6-events, 
which are not considered in the second multiplicative correction. Therefore the remaining 
contribution due to 6-events is taken from the Monte Carlo simulation and subtracted 
from the matrix corrected distribution. This correction is small, the mean of the multi- 
plicity distributions is changed by only ~ 0.6%. 

The central values for the results of the mean multiplicities as a function of the event 
topology are taken to be the arithmetic mean of the respective multiplicity distributions. 
However, in order to estimate the systematic effects of the correction procedure, two 
variations of the correction procedure are also applied: 

• The mean multiplicity is determined from the expectation value of a negative bino- 
mial distribution fitted to every multiplicity distribution. 
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• Instead of correcting the distributions, the means of the uncorrected distributions are 
taken and corrected with a multiphcative factor taken from Monte Carlo simulation. 

The corrected mean multiplicities obtained by using all three correction procedures agree 
well, the small variation is considered in the systematic errors of the analysis. 

This procedure provides fully corrected mean multiplicities of hadronic three-jet events 
with initial udscb-quark - and udsc-quarks for cross-check reasons - as a function of 
the event topology for general and symmetric topologies. Additionally, the multiplicity 
distribution and mean multiplicity of all accepted events regardless of angular restrictions 
is treated in the same way to compare this measurement with previous measurements of 
this quantity. The events entering this overall multiplicity distribution are not required 
to fulfil the cuts on the jet structure given in the lower part of Tab. 2. 



3 Theoretical Predictions 



The ratio of gluon and quark multiplicities r has been subject to theoretical studies 
for some time. The naive expectation that this ratio reflects the ratio of the effective 
colour charges of quarks and gluons given by the colour factor ratio Ca/Cf is subject to 
large corrections. The multiplicity ratio 



99 



N, 



ro[l -ri7o -r27o - ^37o] 



is given as an expansion in the anomalous dimension 70, which can be expressed in terms 
of the strong coupling as as 



7o 



TT 



(2) 



Here Tq denotes the colour factor ratio Ca/Cf, the coefficients of the correction terms 
have been calculated in Leading Order (LO) (ri) [26], Next-to-Leading Order (NLO) 
('"i, ^2) [2] and, using a different approach in which energy conservation is considered, in 
Next-to-Next-to- Next-to-Leading Order (3NL0) (ri, r2 and r^) [27]. As non-perturbative 
effects of the hadronisation process or energy conservation lead to large corrections for 
the ratio r, the suggestion has been made rather to study the ratio of the derivatives of 
the multipficities with respect to the energy scale r^^^ [7] which should be less affected by 
non-perturbative effects. The ratio has been calculated in 3NL0 [3] as 



,(1) ^ d<Ngg> /ds ^ r_ 

d < Nqq > /ds Pi 



(3) 



with 



'2r2ai 
ri 



Po 2 



8C 



+ 



2r2 

oi + ri H J 7o 

ri 



+ 



3r3 



+ ain + Sr2 + — + a2 + aj + rf + 



70^ 



(4) 



and r from Eqn. 1. The coefficients and Oj have been calculated in [3] and are given in 
Tab. 4 and Tab. 5, /3o and Pi are the first to coefficients of the QCD /3-function. 

In another study [15] the ratio r^^^ has been derived in MLLA within the framework 
of the colour dipole model. It is implicitly expressed in the energy evolution of the gluon 
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multiplicity in relation to that of the quark multiplicity: 

dN,,{L' 



dU 



Ca aoCr^ d 

Nqq{L) (5) 



=L+c„ 



Cf V L J dL 



with 



r , / s \ 6 11 3 10 2 3 

L = In — , Q!0 = — , Cg = — , Cg = - , = — TT 



A2y ' ^ 11-2Nf/Ca ' ^ 6 ' ' 2 ' " 27 2 

The constants Cg and Cg are corrections to the phase space available for the quark and 
gluon evolution, Cr determines a MLLA correction which is calculated with an estimated 
uncertainty of 10% and A is the QCD scale parameter. The solution of this differen- 
tial equation implies a constant of integration. Extrapolating the solution of Eqn. 5 to 
small scales, neglecting the constant of integration, would imply that the multiplicity 
in a gg-system would still be significantly larger than in a qq-system. At very small 
scales, however, the hadronic multiplicity of both systems should mainly be determined 
by hadronic phase space and thus should become almost equal [15]: 

NggiLo) ^ NggiLo) = N{Lo) . (6) 

Thus a non-perturbative constant term appears in the solution for the gluon multiplicity 
as expected from the behaviour of the fragmentation functions. In [15] it is suggested to 
determine N[Lq) from data on charmonium or bottomium states. 

The multiplicity of three- jet events is then given by the two alternative formulations 
[16]: ^ 

Nggg = Ngg{Lgg, Ki^,) + -Ngg{KLe) , (7A) 

Nggg = Ngg{L, ^Lu) + ^Nggi^J^^) , (7B) 

henceforth referred to as predictions Eden A and Eden B, with 



and 

The predictions Eqn. 7 A and Eqn. 7B differ in the definition of the gluon contribution 
to the event multiplicity. In Eqn. 7A the ^'^-contribution is determined mainly by the 
invariant mass of the qq-system which is also the relevant scale in an qq-event of the 
same topology with the gluon being replaced by a hard photon. In Eqn. 7B the qq- 
contribution is given by the centre-of-mass energy of the whole event reflecting the phase- 
space available to the g^-pair if no hard gluon had been emitted. 

The expression Ngg{L, n) for the quark contribution to the three-jet multiplicity takes 
into account that the resolution of a gluon jet at a given pt implies restrictions on the 
phase space of the quark system. This restricted multiplicity is linked to the multiplicity 
of an unrestricted qq-system Nqg{L) via [15]: 



Nqg{L, Kent) = Ngq{Kcut + Cq) + {L - Kent " Cq) ^^'^'^^^ ^ 



dU 



(8) 

i'=Kcut+Cq 



7 







02 


0-3 


3 


0.280 


-0.379 


0.209 


4 


0.297 


-0.339 


0.162 


5 


0.314 


-0.301 


0.112 



Table 4: The coefficients aj from 
[3] 





ri 


^2 


rs 


3 


0.185 


0.426 


0.189 


4 


0.191 


0.468 


0.080 


5 


0.198 


0.510 


-0.041 



Table 5: The coefficients from 
[3] 



It is important to note that the unrestricted multiplicity is always a function of only 
one logarithmic energy scale, while the restricted multiplicity demands two logarithmic 
arguments. The equivalent restriction for the gluon contribution to the event multiplicity 
can be neglected, as the evolution scale of the gluon is the p^-like variable k, and therefore 
coincides with the cut-off scale in Eqn. 8. Note also, that the topology dependence of the 
qq-term in Eqn. 7B enters only due to this phase space restriction. 

The explicit energy dependence of the multiplicity of unrestricted ("unbiased") qq 
events has already been calculated in [1] : 



{N{Q')) = a-a':{Q')-exp 



1 + O [V^AQ^) 



with 



1 2Nf 
4^3^ 



1 - 



Ca 



and 



c = 



'327r • Ca 



PI 



(9) 



(10) 



(11) 



In a more recent publication [3] the multiplicity of qq or gg colour singlet systems has 
been calculated in 3NL0 as an expansion of r and 7 = (Ng)' /{Ng), where the prime 
denotes a derivative with respect to the logarithmic energy scale y = log(pO/Qo) with 
being the opening angle of the first splitting and Qq the cut-off of the perturbative 
expansion. The mean multiplicities are then given by 

{Ng{y)) = kg ■ y-'^^^' ■ exp [2C^ + 5g{y)] (12) 

and 

{N,{y)) = - ■ y-"'"'' ■ exp [2CVy + Km . (13) 

kg and denote free normalisations, the quantities Tj and arc calculated in [3] and given 
in Tab. 4 and Tab. 5, C — ^/ANc/Pq and the additional contributions to the exponent 
are given as 



c_ 



2a,C' + § {log(2y) + 2}! + — Lc^ - ^ {log(2y) + 1} 



and 



y 



131 



r2 + 



(14) 



(15) 



Due to the definitions of the evolution variables, which are not symmetric with respect to 
quarks and gluons, the multiplicity of the two-gluon system is given in a slightly higher 
order than the multiplicity of the quark-antiquark system. 
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50 100 50 100 50 100 50 100 



Figure 1: The multiplicity difference between udsc- and wrfscfe-events. The shaded area 
reffects the value for Nq from Eqn. 19, the dashed lines indicate statistical errors only. 

4 Mean multiplicity of hadronic events 

The mean multiplicity of hadronic events at ^/s = mz, without consideration of the 
event topology, is found in this analysis to be 20.963 ±0.004 for udsch- and 20.353 ±0.004 
for Mc?sc-events with statistical errors only. The multiplicity of udsch-evenis is expected 
to be slightly larger than the multiplicity of wdsc-events due to the additional multiplicity 
produced in B decays. These values are in good agreement with previous measurements 
of this quantity [4] . From this measurement the difference in multiplicity between udsch- 
and Mfisc-e vents, which proves to be crucial for this analysis, has been found to be 



udscb—udsc 



udscb 



udsc 



0.610 ±0.002, 



(16) 



where the error of the difference has been determined assuming a correlation of ^1 — Rt ~ 
0.9 between the two measurements. Rb denotes the fraction of hadronic events with initial 
6-quarks, which has been measured on the Z resonance [28]: 



Rb{mz) = 0.21638 ± 0.00066 



(17) 



From previous measurements of the multiplicity of udscb- and 6-events [29,30], Nq can be 
calculated using the relation 



Rh 



Rh 



{nh)b - {rih) 



(18) 



where (n/j)^ and (rih) represent the multiplicity per hemisphere in a 6- or a general 
hadronic event, respectively. The value obtained from the results in [29] is A^^o = 
0.60 ± 0.06, while the value obtained from the results in [30] is A^o = 0.62 ± 0.07, both 
in very good agreement with the value of this analysis. The given statistical errors are 
larger than the statistical error of A^o from this analysis, as in [29] and [30] Nb is measured 



9 




60 80 100 60 80 100 



Figure 2: The mean charged multiphcity of udscb-events with general topologies as a 
function of the opening angle 6i for different values of 9-^ in comparison with the pre- 
dictions of Monte Carlo models. The data points have been plotted in two diagrams 
with alternating ^3-bins, as the neighbouring curves would otherwise overlap. The errors 
shown are statistical only. 

instead of Nudsc, so the determination of Nq is less direct and based on a smaller event 
sample. The systematic errors of Nb are ±0.22 in [29] and ±0.24 in [30]. The average of 
these two values for Nq is 

A^o = 0.61 ± 0.24 (19) 

where also the given systematic errors on A^^ have been considered. 

The additional multiplicity due to initial 6-quarks is in MLLA predicted to be inde- 
pendent of the centre-of-mass energy due to energy conservation and coherence effects in 
the gluon-radiation from heavy quarks [31]. This energy independence has been experi- 
mentally verified (see e.g. [32]). In Fig. 1 the multiplicity difference between udscb- and 
udsc-events Nq is shown as a function of the event topology. The shaded area reflects the 
value given in Eqn. 19, the dashed lines indicate the purely statistical error on A^o- No 
significant tendency can be observed, A'o can be considered as independent of the event 
topology within the angular regions studied. 

The dependence of the multiplicity of M(isc6-events on the event topology is shown in 
Fig. 2 for general and in Fig. 3 for symmetric event topologies. For symmetric topologies 
a nearly logarithmic increase of multiplicity with increasing 9i can be observed. Also 
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Table 6: The multiplicity of udscb-events in dependence of the event topology for general 
topologies. The values for 62 and 6^ represent the centre of the bins, not averages, i.e. due 
to bin overlap it can happen that the given value for 62 is larger than the corresponding 
value for ^3. The first errors are statistical, the second errors are systematic errors due 
to the choice of the cluster algorithm and variations of track and event cuts. 
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Figure 3: The mean charged multiphcity of udscb-events with symmetric topology as 
a function of the opening angle 6i in comparison with the predictions of Monte Carlo 
models The errors shown are statistical only. 

for general topologies the event multiplicity increases with 9i for a fixed value of $3. 
For general event topologies also an increase of multiplicity with decreasing $3 can be 
observed. This dependence is more pronounced than the ^i-dependence as a change in 
61 of 50° results in a multiplicity change of roughly 2 units, while the same change in ^3 
leads to a change in multiplicity of roughly 4 units. The Monte Carlo models, although 
mutually agreeing well, underestimate the multiplicity by ~ 0.4 units. In Fig. 2 it can be 
seen that the deviation is stronger for large values of 63 than for smaller values. 



5 Preparation of the prediction 

Eqn. 7A and Eqn. 7B give the average multiplicity of a three-jet event in terms of the 
restricted gg-multiplicity and the multiplicity of a two-gluon system. Both contributions 
can be derived from the unrestricted gg-multiplicity using Eqn. 8 and Eqn. 5, which has 
been measured in numerous e"^e~ experiments at different centre-of-mass energies. The 
aim of this section is to provide a closed parametrisation of Nqqg using a parametrisation 
of the measured unrestricted multiplicity Nqq{^/s). 

The measurements of Nqq{y/s) used [4] are shown in the central part of Fig. 4. Not 
included are older measurements of the Jade and Pluto collaboration in which the pions 
from the decay — > nn have not been added to the multiplicity [33]. Recent results 
of a re- analysis of the Jade data [34] are, taken into account. In these measurements 
events with initial 6-quarks are included. Due to the decay products of the 6-containing 
mesons these events have an increased multiplicity. As the fraction of events with initial 
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Figure 4: Top: The branching ratio Rh as a function of ^/s Middle: The measurements 
of Nqq [4,5] with the fitted parametrisation Eqn. 22. Note that the 4th order polynomial 
P4 fit does not account for the additional multiphcity due to 6-events. Bottom: The 
relative deviation of the fitted polynomial P4 from the parametrisation Eqn. 9 
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j\^udscb 


^1 


ATUdSCb 


22° 


20.77 ±0.03 


±0.24 


72° 


27.11 


±0.13 ±0.10 


27° 


21.91 ±0.04 


±0.18 


77° 


27.40 


±0.14 ±0.08 


32° 


22.70 ±0.05 


±0.13 


82° 


27.71 


±0.15 ±0.11 


OTO 

37 


23.48 ±0.06 


±0.10 


87 


28.05 


±0.16 ±0.06 


42° 


24.11 ±0.07 


±0.09 


92° 


28.21 


±0.17 ±0.13 


47° 


24.63 ±0.08 


±0.09 


97° 


28.44 


±0.18 ±0.10 


52° 


25.30 ±0.09 


±0.06 


102° 


28.89 


±0.19 ±0.14 


57° 


25.81 ±0.10 


±0.08 


107° 


28.76 


±0.19 ±0.07 


62° 


26.32 ±0.11 


±0.10 


112° 


29.14 


±0.20 ±0.07 


67° 


26.67 ±0.12 


±0.05 


117° 


28.68 


±0.22 ±0.14 



Table 7: The multiplicity of udscb-events in dependence of the event topology for sym- 
metric topologies. The first errors are statistical, the second errors are systematic errors 
due to the choice of the cluster algorithm and variations of track and event cuts. 



6-quarks varies with the centre-of-mass energy, this additional multiplicity introduces an 
additional energy dependence of the multiplicity which is not due to strong interactions 
and which has to be corrected for. The upper plot of Fig. 4 shows the branching ratio 

^ a{e+e- ^ bb) 
a{e+e~ hadrons) 

as a function of the ccntrc-of-mass energy. The curve has been obtained using the LUX- 
TOT routine of the Jetset/Pythia package [24]. The average additional multiphcity 
due to 6-events can be expressed in terms of Nq: 

Sb-udsc = 2 • ({nh)b - {nh)udsc) = D ^ ■ (21) 

With the values from Eqn. 17 and Eqn. 19, the multiplicity difference can be determined 
giving Sb-udsc = 2.83 ± 0.23. 

The gg-multiplicities are then fitted with the function 

Nudscb{^/s) = Nudsc{^/s) + Sb-udsc ■ Rb{^/s) (22) 

where Eqn. 9 is used to parametrise NudsdVs) ■ The description of the data is very good 
with a per degree of freedom of 33/48. The values found for the fit parameters are 

a = 0.10252 ± 0.0025 

A = (0.243 ± 0.012) GeV . 

The fitted function is indicated in Fig. 4 as a solid fine. To simphfy the integration of 
Nudsc(\/s) required due to Eqn. 5, a polynomial of order four in L = log(s/A^) is fitted 
to the parametrisation of Nudsci^/s)- The fit is performed with A = 0.25 GeV, and the 
values obtained for the coefficients are given in Tab. 8. The polynomial is indicated as a 
dashed fine in the middle plot of Fig. 4. The difference between the dashed polynomial 
and the solid line indication of Eqn. 22 is due to the omission of the Sb-udsc • RbiV^) 
term. In the lower part of Fig. 4 the relative deviation of the polynomial from the fitted 
Nudsc{^/s) is shown. The deviations are smaller than 2 ■ 10~^ over most of the fitted 
region and nowhere exceed 4 • lO"'', indicating an excellent description of Nudsc{y/s) by 
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Po 0.18981 
pi 0.40950 
P2 0.57358 
P3 0.10349 
P4 0.28640 



10-1 
10-2 
10-3 



Table 8: The coefficients of the polynomial fitted to N^h^^{-\/s) 



the polynomial. The value of the polynomial at y/s — mz is 20.2561. This value is in 

good agreement with the mean multiplicity of wrfsc-events of 20.353 ± 0.004 measured in 
this analysis. The small difference is taken into account when considering the systematic 
errors of this analysis. 

Inserting the polynomial into Eqn. 5 and Eqn. 8 leads to a closed form of the predic- 
tions of Eqn. 7A and Eqn. 7B. To fix the constant of integration, left free in Eqn. 5, a 
measurement of A^^^^ from x' decays by the Cleo collaboration is used [35]. The value 
of Ngg{Qm?>2 GeV) = 9.339 ± 0.090 ± 0.045 using A = 250 MeV results in Lq = 5.86, 
which corresponds to a centre-of-mass energy of 4.68 GeV at which quark and gluon 
multiplicities are equal. 

Assuming massless jet kinematics, the scale variables Km, i^hc and Lqq can directly 
be expressed only in terms of the inter-jet angles di (and the constant centre-of-mass 
energy). This yields Nqqg as a function of the inter-jet angles under the assumption that 
a certain jet is the gluon jet. As no explicit identification of the gluon jet is made, Nqqg 
is calculated for all three possible gluon jet hypotheses and the weighted mean of these 
values is taken, where each hypothesis is weighted with the corresponding value of the 
three-jet matrix-element which can be calculated from the inter-jet angles assuming again 
massless jet kinematics. 

However, it is a known feature of the hadronisation process that due to coherent pro- 
duction of inter-jet particles close-by jets are pulled even closer together. The calculations 
leading to Eqn. 7A and Eqn. 7B refer to the partonic structure of an event. While for 
the change in multiplicity due to the hadronisation process in accordance with the LPHD 
hypothesis an overall normalisation constant can be found, for the change in the event 
topology, i.e. the inter-jet angles, this is not the case. Therefore, a topology dependent 
hadronisation correction has to be applied. The effect can be most easily understood 
for symmetric event topologies. During the hadronisation process jets two and three are 
pulled closer together resulting in a smaller opening angle 6i at the hadronic than at the 
partonic level of the event. As the multiplicity increases with 9i, the reduced 9i found at 
the hadronic level thus results in an underestimation of the multiplicity compared to the 
prediction which is based on the value of 9i at parton level. 

This effect is corrected for by a topology dependent correction factor which is applied 
to the prediction. In Fig. 5 the correction factors obtained for symmetric event topologies 
arc shown. The correction factors are calculated by dividing the mean hadron multiplicity 
obtained for a certain 9i bin (with 9i measured at the hadronic level of the simulated 
event) by the mean hadron multiplicity obtained with 9i measured at the partonic level 
of the event. The correction factors shown are taken from the Ariadne Monte Carlo 
simulation for each of the four cluster algorithms used. The correction factor is around 
unity for large opening angles and increases for small 9i. This follows the expectation 
that jets which are close together are pulled even closer together than jets with a larger 
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Figure 5: The hadronisation correction factor for symmetric topologies from the Ariadne 
Monte Carlo. Plotted is the ratio of multiplicities with 9i determined at the hadronic 
level and the multiplicities with Oi determined at the partonic level. 
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Table 9: The parameters of the hadronisation correction factors for symmetric event 
topologies 



angle between them. In order to get a smooth correction the function 

c = (23) 

is fitted to the correction factors. The correction is described very well by the fitted 
function, the corrections obtained for the four studied cluster algorithms are shown in 
Fig. 5. The fitted parameters obtained using the Monte Carlo generators Ariadne, 
Herwig and Jetset are given in Tab. 9 for the four cluster algorithms. The same 
procedure is applied also to events with general topologies. In Fig. 6 the correction 
factors obtained with Ariadne for the four used cluster algorithms are shown for several 
values of O^, as a function of 6i. In order to describe these corrections smoothly the 
two-dimensional function 

c = (oi + 02^3 + 05^3) + (03 + 04^3 + Qe^l) • 0\ (24) 

is fitted to the correction factors. The description of the correction by the fit is reason- 
able, and the parameter values obtained are listed in Tab. 10. In general, the correction 
becomes smaller with increasing opening angle ^1, as expected for this hadronisation 
effect. 
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Figure 6: The hadronisation correction factors for general topologies. Plotted is the 
ratio of multiplicities with the inter- jet angles determined at the hadronic level over the 
multiplicities with the inter-jet angles determined at the partonic level. 
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Table 10: The parameters for the hadronisation correction factors for general topologies 
for the angular ordered Durham (a.o.D.), Cambridge (Camb.), Durham (Durh.) and 
Lucius (Lucl.) algorithm 

6 The fit of CaICf 

The closed form of Nqqg as a function of the event topology obtained from Eqn. 7A 
(Eden A) and Eqn. 7B (Eden B) with the procedure described in Sect. 5 can now be 
compared with the measured three-jet event multiplicities. In Fig. 7 the measured multi- 
plicities of udsc- (open markers) and tidscfe-events (solid markers) are shown for symmetric 
and general topologies. The two solid lines indicate the respective predictions of Eqn. 7A 
where for rtdsc6-events the constant A^o of Eqn. 19 has been added to the prediction. The 
dashed lines represent the respective predictions of Eqn. 7B. It can be observed that pre- 
diction Eqn. 7A describes the multiplicity of symmetric events in an excellent way, while 
Eqn. 7B overestimates the multiplicity by ~ 0.6, with the result that the prediction for 
wdsc-events almost coincides with the multiplicities measured in rtdsc6-events. Moreover, 
the slope of Eqn. 7B seems to be larger than the slope of the measured multiplicities. 
For general topologies an overall good description of the event multiplicities by Eqn. 7A 
can be seen. Only for large 6*3 and then especially for large 61, i.e. in topologies where jet 
3 is not strongly pronounced, larger deviations occur. Again, Eqn. 7B overestimates the 
multiplicities significantly and shows a larger slope with respect to 61. 

In order to determine the colour factor ratio Ca/Cp from the event multiplicities, the 
predictions Eqn. 7A and Eqn. 7B are fitted to the data. To avoid systematic uncertainties 
entering through the use of a 6-tagging procedure, M(isc6-events are used to obtain the 
central result, while the results obtained for itdsc-events are considered when estimating 
systematic uncertainties. Ca/Cf enters the prediction only via the derivative of the 
gluon multiplicity in Eqn. 5, so this parameter is expected to be sensitive to the change 
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Figure 7: The mean multiplicity of udscb- and udsc-events with symmetric and general 
topologies in comparison with the predictions Eqn. 7A and Eqn. 7B prepared as described 
in Sect. 5. The error bars indicate statistical errors only. 
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Figure 8: The fit parameters x^/^df, Ca/Cf and Nq as function of the fit range. The 
shaded area indicates the expected value of A^^o with statistical and systematical uncer- 
tainties, the dotted lines indicate statistical errors only. The dashed lines indicate 1.0 in 
the top, Ca/Cf = 2.25 in the middle and the A^o value from Eqn. 19 in the bottom plots. 
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udscb 
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2.257±0.019 


0.62±0.03 


2.127±0.034 


0.12±0.03 


gen. 


2.261±0.014 


0.60±0.03 


2.028±0.050 
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Table 11: The results of the parameter fit for Ca/Cf and A^o 



of multiplicity rather than to its absolute value. Therefore, the additional constant A^^o 
is allowed to vary freely in the fit (also when fitting udsc-events where A'o is expected to 
be zero) in order to avoid an influence of the absolute value on the flt result for Ca/Cf- 

The flt range used is determined by the value of X^/^df obtained with the flt. While 
for symmetric event topologies it is straightforward to determine the data points used 
in the fit by demanding a minimum opening angle 6i, the situation is not as clear for 
the general event topologies, where two variables are needed to describe a topology. A 
reasonable order criterium is a 'three-jet likeliness' of a topology, i.e. the emphasis which 
is put on the third jet. Therefore the pt of the third jet is a reasonable ordering variable. 
It turned out that Pt,Lu is the best suited of the used p^-like variables, because the fits 
converge fastest, if the fit range is determined by Pt.Lu- The fit range is then defined by 
demanding a minimum pt^Lu-value for the event topologies used. 

In the upper plots of Fig. 8 the values obtained for x^/A^/ when fitting udscb-events 
with the predictions of Eqn. 7A and Eqn. 7B are shown for both general and symmetric 
topologies. The values of /N^f obtained for symmetric events are always smaller than 
1.0, regardless of the chosen minimal ^i-value except for extremely small fit-ranges. The 
curves show a shallow minimum around ~ 30° which is chosen to determine the 
fit-range. For general topologies, and a small minimal Pf,Lu; the X^/^df values are large 
but decrease fast with a more restricted fit range. The /Ndf curve shows a step-like 
structure, where every step corresponds to a completely excluded bin in 6*3. From a 
minimum pt^hn value of ~ 22 GeV on a plateau in the vicinity of x^/^df — 1 can be 
observed and a minimum pt,hu of 25 GeV is chosen. 

Fitting the data within these ranges results in the parameter values given in Tab. 11. 
As can be seen, the prediction of Eqn. 7A results in compatible values for Ca/Cf when 
fitting udsc- and udscb-events of symmetric or general topologies. Both values for A^q 
obtained by fitting udscb-events agree with the expectation of Aq = 0.61 as well as both 
Ao values for itdsc-events agree with the expectation of Aq = 0. The results obtained for 
Ca/Cf and A'q with Eqn. 7B in symmetric topologies are not in good agreement with 
the values obtained in general topologies. Also the values for Aq are not in agreement 
with the respective expectations. In the lower plots of Fig. 8 the fit results for Ca/Cf 
and A^o are shown for udscb-events as a function of the fit range. The parameters are 
correlated, but as this plot is only to show the influence of the flt range, when plotting 
the dependence of one parameter the other parameter is fixed to the respective value 
of Tab. 11. While Eqn. 7A leads to results which are independent of the fit range, the 
results of Eqn. 7B exhibit a more pronounced dependence on the fit range. Again, the 
discrepancy of the results obtained with Eqn. 7B for symmetric and general topologies 
can be observed, as well as the disagreement of the obtained values for A'o with the 
expectation indicated by the shaded area while fixing Aq to the expected values leads to 
an unacceptably high /Ndf- 
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variable 


tight 


loose 


P 


> 0.5 GcV 


> 0.3 GeV 


^polar 


30° - 150° 


20° - 160° 


^xy 


< 2 cm 


< 8 cm 




< 8 cm 


< 12 cm 




< 80% 


< 120% 



Table 12: Variations of the track cuts 
given in Tab. 1 for systematic studies 



variable 


tight 


loose 


general events 


pihemisph. 

charged 


> 0.05- 


> 0.02 ■ ^ 


ptotal 
charged 


> 0.18- ^ 


> 0.06 • ^ 


-^charged 


> 6 


> 4 


^sphericity 


40° - 140° 


30° - 150° 


Pmax 


20 GeV 


70 GeV 


three-jet events 




> 357.5° 


> 350° 


-E'visible/jet 


> 8 GeV 


> 3 GeV 




40° - 140° 


30° - 150° 



Table 13: Variation of the event cuts 
given in Tab. 2 for systematic studies 



Since the prediction Eqn. 7B, which depicts a rather extreme scenario of the phase- 
space assignment to quark and gluon jets, fails to give consistent results for symmetric 
and general topologies and because Eqn. 7B overestimates the event multiplicity in gen- 
eral which results in unphysically low values of A^o, this formulation of the prediction is 
discarded in favour of Eqn. 7A for the central result. The central results of the fits are 

2.257 ±0.019stat. (25) 

2.261 ± 0.014,tat. (26) 

for general event topologies. 

Several sources of systematic uncertainties have been considered. The systematic error 
on Ca/Cf is estimated by fitting prediction Eqn. 7A with the variations to the analysis 
listed below. As the parameters Ca/Cp and A^o ^-re correlated and this correlation should 
not be reflected by the estimated systematic uncertainties, A'o is flxed at 0.61. The values 
obtained for Ca/Cf with the variations of the analysis are then compared to the result 
obtained by the standard analysis but with A'o fixed at the above value. 

The studied experimental sources of uncertainty are: 

• Variations of the track cuts hsted in Tab. 1 within the ranges given in Tab. 12. The 
minimal track length is not lowered in order to maintain a contribution of the TPG 
detector to the track reconstruction; 

• Variations of the cuts on the event- and jet-structure listed in Tab. 2 within the 
ranges given in Tab. 13 ; 

• The hadronisation corrections are calculated from Jetset instead of Ariadne sim- 
ulation. Alternatively 30% of the hadronisation correction are regarded as this 
correction's uncertainty; 

• Alternatively to the arithmetic mean of the multiplicity distributions, the mean 
multiplicities are determined by parameter fits of negative binomials to the distribu- 
tions. As a further alternative, the matrix correction is replaced by a multiplicative 
correction to the mean multiplicities; 

• The udsc sample is used to estimate the uncertainty of Ca/Cf due to the 6-quark 
multiplicity; 
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0.4% 
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0.1% 
0.5% 
0.9% 
0.8% 



variation of A 
variation of Cr 
variation of Lq 
clustering algorithm 



1.0% 
2.0% 
0.0% 
2.2% 



2.5% 



4.0% 



3.1% 



0.4% 
0.6% 



0.4% 
0.2% 
0.8% 
0.0% 
0.5% 
0.7% 
0.7% 



1.0% 
2.1% 
0.0% 
1.8% 



1.6% 



2.9% 



3.3% 



Table 14: The relative systematic uncertainties of the measurement of Ca/Cf 



• The lower limit of the fit range is varied between 6i > 25° and 6i > 35° for symmetric 
events and between pt^hu > 22.5 GeV and pt^hu > 28 GcV for general topologies; 

• The difference between the parametrisation of Ngg at y/s = rriz and the event mul- 
tiplicity measured in this analysis is 0.101 units, i.e. 0.5%. In the worst case this 
reflects an uncertainty on the normalisation of A^^^^ which would directly be reflected 
in Ca/Cf- Therefore this deviation is considered conservatively as the relative un- 
certainty on Ca/Cf] 

• The data set used to fit the parametrisation of Nqq is varied, by excluding the 
results of the ToPAZ, Jade and Mark-I collaborations or the results of the Amy 

collaboration respectively; 

• The correction due to the additional multiplicity in 6-events is varied within the 
errors of 5b^udsc when fitting the parametrisation of Nqq. 

Additionally, the uncertainty due to the choice of the clustering algorithm has been 
estimated by comparing the result obtained with the angular ordered Durham algorithm 
with the results obtained using the Durham, Cambridge and Lucius algorithms. Moreover 
the following uncertainties of the theoretical prediction have been studied: 

• The scale variable A has been varied between 200 MeV and 300 MeV; 

• The constant Cr has been varied within the given theoretical uncertainty of ~ 10%; 

• The measurement of Ngg by Cleo which has been used to determine Lq has been 
varied within the given errors. As Lq is a constant of integration it is supposed to 
affect mainly Nq and not Ca/Cf- Therefore Nq is varied freely when estimating the 
effect of variations of Lq. Indeed no dependence of Ca/Cf on the variation of Lq is 
observed. 

The relative uncertainties caused by the different sources are given in Tab. 14. Alter- 
natively to estimating the systematic uncertainty inherent to the prediction Eqn. 7A by 
varying the parameters as discussed above, it is suggested in [36] to study the difference 
between the prediction Eqn. 5 based upon the colour-dipole model with predictions based 
upon parton-shower evolution, as e. g. Eqn. 3. Since this measurement of Ca/Cf is sen- 
sitive rather to the energy dependence of the multiplicities than to the absolute value 
of multiplicity, the uncertainty of Eqn. 5 is relevant here. Both predictions, Eqn. 3 and 
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/^Le [GeV] 



Ngg ±(stat.)±(sys. 



«Lc [GeV] 



Ngg ±(stat.)±(sysO 



13.36 
15.57 
17.82 
20.13 
22.50 
24.92 
27.45 
30.01 
32.61 



11.21± 0.14 
12.41± 0.12 
13.43± 0.11 
14.32± 0.11 
15.56± 0.11 
16.53± 0.11 
17.57± 0.12 
18.32± 0.13 
19.27± 0.14 



± 0.14 
± 0.08 
± 0.09 
± 0.08 
± 0.08 
± 0.13 
± 0.16 
d= 0.10 
± 0.15 



35.26 
38.01 
40.69 
43.37 
45.87 
48.19 
50.19 
51.69 
52.50 



19.96± 0.14 
20.73± 0.16 
21.56± 0.16 
22.09± 0.17 
22.72± 0.18 
23.81± 0.19 
23.71± 0.19 
24.60± 0.20 
23.75± 0.22 



± 0.11 
± 0.13 
± 0.15 
± 0.18 
± 0.18 
± 0.21 
± 0.11 
± 0.15 
± 0.22 



Table 15: Ngg extracted from three-jet events with symmetric topologies 



Eqn. 5 are shown in Fig. 13. The difference between the two predictions, evaluated at 
the centre of the data around ~ 40GeV is ~ 3.7%. This corresponds to an one-sigma 

error-estimate of 2.1% which is completely compatible with the theoretical error given in 
Tab. 14. The result obtained with symmetric topologies is therefore 

^ = 2.257 ± 0.019stat. ± 0.056exp. ± 0.070theo. (27) 

and 

§i = 2.261 ± 0.014,tat. ± 0.036exp. ± 0.066theo. (28) 

with general topologies. Both results are strongly correlated due to a large fraction of 
common events, so an average cannot be made here. Instead, the more precise result 
obtained with general topologies is taken as the central result. This result is the most 
precise measurement of the colour factor ratio so far with an overall relative uncertainty 
of 3.4%. 

In Fig. 9 the colour factor ratios Ca/Cf and Tr/Cf are mapped for several symmetry 
groups where = NpTp with A^^ as the number of active quark flavours and Tp = 1/2 
as the normalisation of the SU(3) representation. In this plot the present analysis result 
is indicated by the shaded vertical band. The shaded diagonal band is the result of a 
measurement of the QCD /3-function correlating Tr/Cf with Ca/Cf [37]. The dashed 
ellipses represent measurements of the four-jet cross section as function of the inter-jet 
angles [38]. The measurements are combined by adding their x^-functions. The solid 
contour represents the A^^ = 2.4 limit, corresponding to a confidence level of 68%. It is 
in excellent agreement with the QCD expectation of SU(3). 



7 The extraction of Ngg 

Instead of fitting Ca/Cf with the measured three-jet event multiplicities, the predic- 
tions of Eqn. 7A and Eqn. 7B can be used to extract the multiphcities of unrestricted 

two-gluon colour-singlet systems from the three-jet event multiplicity. As the direct ex- 
perimental access to this quantity is severely limited, this proves to be an interesting 
option. The predictions Eqn. 7A and Eqn. 7B are solved for the gluonic contribution 
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1 2 



c /c 

Figure 9: The plot shows the Casimir eigenvalues expected for the special orthogonal 
(SO(N), straight line on the left) and the special unitary (SU(N), curved hne on the 
right) groups as well as for the symplectic groups (SP(2N), straight line on the right). 
Eg, E7, Eg, G2 and F are the five exceptional Lie groups. The shaded bands and the 
dashed ellipses indicate the result of this and other [37,38] experimental analyses, which 
are combined in the solid ellipse. 

yielding 

Nggini^e) = 2 ■ (iV,,-3(t92, ^3) - Ng^iL^g, ^Lu)) (29A) 

NggiKi^^) = 2 ■ (^N,^gi^2, ^3) - N,^{L, /€Lu)) . (29B) 

As discussed already in the previous section, Eqn. 7B fails to describe the data so only 
Eqn. 29A, which is based on Eqn. 7A, is used to extract the gluon multiplicity. The 
hadronisation correction discussed in Sect. 5 is taken into account by dividing Nqqg by 
the appropriate correction factor. Nqq is calculated for all three gluon-jet hypotheses and 
averaged with the appropriate weights as discussed above before subtracting it from the 
three-jet multiplicity. Analogously, the value of klg, at which Ngg is being determined, is 
the weighted average of the ^Le-values for the three possible gluon-jet hypotheses. 

According to the procedure of Sect. 6 the angular ordered Durham algorithm is applied 
to obtain the central results. Also, only the topologies which entered the fit in Sect. 6 
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Table 16: A^^^ extracted from three-jet events with general topologies 



are considered in order to provide a good description of the three-jet multiplicity by the 
prediction. As the additional multiphcity due to 6-quark decays has been found to be 
topology independent, udscb-events are used in order to avoid the uncertainties inherent 
in the 6-tagging procedure and Nq is fitted to the data yielding Nq = 0.628 ± 0.017 for 
general and A^"o = 0.628±0.024 for symmetric topologies with the x^/^df values of 80.4/74 
and 7.3/17 respectively. The differences between the prediction and the measured three- 
jet event multiplicity are shown in Fig. 10. The agreement between data and prediction 
is good, especially for symmetric topologies. Anti-6-taggcd w,(i,sc-cvcnts, with Nq fixed at 
0, are taken into account when estimating the systematic uncertainties of Ngg. 

The gluon multiplicities obtained are given in Tab. 15 and Tab. 16. They are shown 
in Fig. 11 as a function oi pt^Le which has been identified with the centre-of-mass energy 
of the unrestricted two-gluon colour-singlet system. Solid circles indicate the results from 
symmetric, solid squares from general event topology. For clarity the multiplicities from 
general topologies have been rcbinncd in pt in the upper plot. The unrcbinncd results 
are shown in the lower plot of Fig. 11. Results of both topology classes are found to 
be in good agreement. The good agreement between prediction and data in symmetric 
topologies also for small opening angles allows gluon multiplicities at small values of /«Le 
to be obtained resulting in a larger kinematic range covered by the gluon multiplicities 
from symmetric than from general topologies. 
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Figure 10: Difference between tlie measured tfiree-jet event multiplicity and the prediction 
Eqn. 7A. 

The solid line in Fig. 11 represents the prediction for the gluon multiplicity which 
has been described in Sect. 5. The agreement between this prediction and the extracted 
gluon multiplicities reflects the good description of the three-jet event multiplicities by 
Eqn. 7A. The measurement of Ngg of the Cleo collaboration at ~ 10 GeV shown in 
Fig. 11 has been used to fix the constant of integration when setting up the prediction 
for Nggj therefore it agrees by definition with the prediction. However, this vahic docs 
not enter in the determination of Ngg. The measurement at ~ 5 GeV is represented with 
statistical errors only. The Opal measurement of Ngg around 80 GeV is based on the 
measurement of gluon jets recoiling against two identified 6-quark jets [10]. The Opal 
measurements at lower energies used identified gluon jets from three-jet events and an 
effective jet energy scale [12]. Also included is a measurement of the Topaz collaboration 
obtained from fully symmetric events [6]. The measurement [18] is not included in Fig. 11, 
as it is based on Eqn. 7B which has been found not to describe the data satisfactorily. 
The overall agreement between the several measurements is good. 

After the preliminary presentation of the analysis of symmetric three-jet events in [17], 
a similar analysis has been pubhshed [18], where different conclusions have been reached. 
Especially the prediction Eqn. 7B has been preferred over Eqn. 7A due to the observed 
extrapolation behaviour of the extracted gluon multiplicities. However, only symmetric 
topologies have been studied in [18], while the main arguments of this analysis to disfavour 
Eqn. 7B compared to Eqn. 7A are the inconsistent results for symmetric and general 
topologies obtained when using Eqn. 7B. Moreover, the extrapolation behaviour of the 
gluon multiplicities extracted with Eqn. 7A to direct measurements at lower energies is 
satisfying, when a hadronisation correction is applied. Note that in [18] no hadronisation 
correction has been used and the extrapolation behaviour of the gluon multiplicities 
extracted using Eqn. 7A has been found to be unsatisfying, leading to the preference of 
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Figure 11: Gluon and quark multiplicities. For clarity, gluon multiplicities from general 
topologies are rebinned in the upper plot. In the bottom plot these Multiplicities are 
shown for each topology bin. 
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Figure 12: Measurements of r as defined in Eqn. 1 in comparison with theoretical pre- 
dictions. The shaded bands indicate the effects of a variation of Np between 3 and 5. 



Eqn. 7B over Eqn. 7A in [18]. A more extensive attempt to compare both analyses can 
be found in [39]. 

In addition to the gluon multiplicity, in Fig. 11 the unrestricted gg-multiplicity mea- 
sured by several e~^e~ experiments which have been used to parametrise the prediction 
are shown. The measurements have been corrected for the varying 6-contributions, so the 
multiplicity of udsc-events is shown. The line indicates the fitted parametrisation from 
Eqn. 9. The higher multiplicity and the larger slope of the gluon multiplicity with respect 
to energy in comparison with the qq multiplicity can be clearly observed illustrating the 
observation of the larger colour charge of gluons. 

In order to obtain the ratio r of the gluon and quark multiplicity, the extracted gluon 
multiplicities are divided by the parametrisation of N^g evaluated at the respective energy. 
The obtained values for r are shown in Fig. 12. Open dots represent the results from 
general, solid dots from symmetric topologies. The curves indicate the LO, NLO and 
3NL0 predictions of r discussed in Sect. 3 with the shaded area representing the effects 
of a variation of Np between 3 and 5. Additionally, the result of a numeric calculation [40] 
is indicated. The predictions clearly overestimate the measured ratio. The overestimate 
is larger the lower the order of the calculation is. The difference between the predictions 
is also rather large. The reason for the poor description of the measurement by these 
predictions are non-perturbative effects to which the absolute value of multiplicity is 
especially sensitive. Only the line derived from the prediction of Eqn. 5 describes the data, 
but in the determination of this curve experimental input, especially the measurement of 
the gluon multiplicity at a/s = 10 GeV which fixes the constant of integration, enters, so 
this line includes a non-perturbative correction. 

The ratio of the derivatives of gluon and quark multiplicity with respect to the energy, 
r^^\ is expected to be less sensitive to such non-perturbative effects, r^^^ is obtained from 
the measured gluon multiplicities by taking the linear slope of mutually exclusive pairs 
of measured multiplicities which is then divided by the derivative of the parametrisation 
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Figure 13: Measurements of r^^^ as defined in Eqn. 3 in comparison with theoretical 
predictions. The shaded bands indicate the effects of a variation of Np between 3 and 5. 





k 


Qo[GeV] 






0.243 ±0.011 
0.1980 ±0.0005 


0.454 ±0.051 


0.75 
8.5 




0.343 ±0.012 


1.317 ±0.084 


1.1 



Table 17: Parameter values of the 3NL0-fits of N„iy and N, 



of Ngg taken at the respective energy value. The obtained values for r'-^^ are shown 
in Fig. 13. Only symmetric topologies have been considered here, as they cover a larger 
kinematic range and are more evenly distributed in energy. However, the results obtained 
from general topologies for r^^^ are in full agreement with the results shown in Fig. 13. 
The measured r^^-* are constant within their errors, the weighted mean is: 

rW = 1.75 ±0.07 . (30) 

r^^^ is therefore significantly higher than r due to the weaker dependence on non- 
perturbative effects, which reduce the ratio from the asymptotic value of Ca/Cf- The 
predictions Eqn. 3 and Eqn. 5 are indicated in Fig. 13 as curves. A reasonable agreement 
between the perturbative calculation and data can be observed here. Note that Eqn. 5 is 
a purely perturbative prediction for the slopes. Also the difference between both predic- 
tions is smaller for r^^^ than for r. Both observations confirm that the energy dependence 
of multiplicity is an observable superior to the absolute value of multiplicity. 

The scale-independent values obtained for r^^^ already indicate that there is no direct 
sensitivity to the ratio r^^^ of the second derivatives of the multiplicities with respect 
to energy. However, the parametrisations of Eqn. 12 and Eqn. 13 can be fitted to the 
data and derivatives can then be obtained analytically. This procedure has already been 
applied in [18] in order to obtain r^'^\ As already observed in [3], the use of a common 
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Figure 14: The top plots show the 3NL0 predictions Eqn. 12 and Eqn. 13 fitted to the 
measured muhiphcities. The lower plot shows from bottom to top r, r^^^ and r^"^^ which 
are derived from the fitted parametrisations together with the measured points for r and 
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normalisation k for quark and gluon multiplicity leads to a significant underestimation of 
the quark multiplicity and a wrong energy dependence. Therefore, the normalisations kg 
and kg are allowed to vary independently. As the quark multiplicities are measured more 
directly, Qo is fixed by the fit of Nqq. This leaves kg as the only free parameter in Ngg. 
The values obtained for the parameters are given in Tab. 17, the fitted paramctrisations 
are shown in the upper two plots of Fig. 14 as a solid line for Ngg and as a dashed line 
for Ngg. The fit of Ngg describes the data well, the value obtained for is small and 
Qo = 454 MeV is of of the same order of magnitude as Aqcd- However, the energy 
dependence of the gluon multiplicity is not described well with this value of Qq, the 
obtained is quite large. Therefore, Qo is also allowed to vary freely in the fit of Ngg 
resulting in two different effective energy scales for quarks and gluons. This apparently 
implausible finding can be motivated by the fact that due to the choice of evolution 
parameters leading to Eqn. 12 and Eqn. 13 both quantities are predicted not in exactly 
the same order leading to different effective scales. The variation of Qq leads to a good 
description also of the gluon multiplicities by the prediction as indicated by the solid line 
in the upper right plot of Fig. 14, the parameter values obtained for Ngg are given in the 
last line of Tab. 17. 

The derivatives of these fitted paramctrisations can now be calculated analytically to 
any order and the ratios of these derivatives can be built. In the lower plot of Fig. 14 from 
bottom to top the ratios r, r^^^ and r^^^ are shown in comparison with the measurements 
of r and r^^\ The agreement between the measurement and the derived ratios is good, 
especially the energy dependence of r is well reproduced. The shaded area indicates the 
result obtained for r*^^^ within the uncertainty given by the errors of the fit parameters. 
Interpreting this curve as a measurement of r*^^-* the value obtained e.g. at 30 GeV would 
be found as 

r^^^ (30 GeV) = 1.97 ± 0.09 . (31) 

This value is in good agreement with the value obtained by a similar procedure in [18]. 
However, theoretical assumptions about the energy dependence of quark and gluon mul- 
tiplicity are the basis of the paramctrisations used, so that a theoretical bias cannot be 
avoided in obtaining this quantity. 



8 Conclusions 

The multiplicity of hadronic thrcc-jct events has been measured as a function of vari- 
ables depending on the event topology. A MLLA-prcdiction of this quantity has been 
fitted to the data yielding a measurement of the colour factor ratio 

^ = 2.261 ± 0.014,tat. ± 0.036exp. ± 0.066theo. 

Up 

With an overall relative uncertainty of 3.4% this is the most precise measurement of this 
quantity so far. It has been shown that the formulation Eqn. 7B of the prediction should 
be rejected in favour of the formulation Eqn. 7A. Using the MLLA-prediction to subtract 
the quark-contribution from the three-jet event multiplicity, the multiplicity of two-gluon 
colour-singlet systems has been extracted over a wide range of the effective energy scale. 
Comparing these multiplicities with the known multiplicities of quark- antiquark colour- 
singlet systems, the ratios r, r^^^ and r^^^ have been studied. It has been found that r^^^ is 
an observable superior to r as non-perturbative effects affect the energy development of 
multiplicity less than the absolute value of multiplicity. However, the data do not show 
any significant sensitivity to the ratio r^'^\ 
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